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ABSTRACT 

An  equation  is  derived  that  approximately  describes  pressure  re¬ 
sponse  in  a  long  constant-diameter  tube  which  is  closed  at  one  of  its  ex¬ 
tremities.  Numerical  solutions  are  presented  for  several  types  of  forc¬ 
ing  fiinctions. 
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SUMMARY 


An  analysis  is  made  of  pressure  response  in  a  system  that  consists  of  a  tube  open  at  one  extremity  and 
terminating  at  the  other  end  in  a  reservoir.  The  analysis  concerns  the  continuum  regime  and  pertains  to  fully 
developed  laminar  flow  of  a  compressible  fluid  with  isothermal  state  changes.  An  equation  is  derived  that  ap¬ 
proximately  describes  pressure  response  in  the  system  when  the  tube  is  long  and  the  reservoir  volume  is  small. 
The  response  equation  is  made  dimensionless,  and  numerical  solutions  are  presented  for  three  types  of 
boundary  conditions:  step  functions,  ramp  functions,  and  the  functions  dictated  by  constant-velocity  flight  in  an 
isothermal  atmosphere.  The  results  of  this  analysis  may  be  used  to  determine  response  in  pressure -sensing 
systems  during  the  transient  conditions  of  flight  in  real  atmospheres,  and  in  systems  subjected  to  impulse  pres¬ 
sure  functions. 
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PRESSURE  RESPONSE  IN  A  LONG  TUBE 
WHICH  IS  CLOSED  AT  ONE  END 


Introduction 

It  is  often  necessary  to  measure  a  thermodynamic  property  such  as  pressure  or  density  during  transient 
conditions.  A  typical  sensing  system  may  consist  of  a  tube  which  is  joined  at  one  extremity  to  an  orifice  and  is 
connected  at  the  other  end  to  a  reservoir  that  has  an  element  suitable  for  measurement  of  the  thermodynamic 
property.  If  the  tube  is  long  or  the  reservoir  volume  is  large,  the  measured  value  may  be  significantly  differ¬ 
ent  from  the  value  applied  at  the  orifice,  particularly  when  the  pressure  is  low.  To  interpret  that  type  of 
measurement  or  to  design  sensing  systems,  the  response  characteristics  of  the  system  must  be  established. 

Vaughn  (see  Bibliography)  derives  an  equation  that  describes  pressure  response  in  a  reservoir  with  a 
tube  of  negligible  length.  Vaughn's  equation,  which  is  nonlinear,  permits  a  more  accurate  estimate  of  pres¬ 
sure  response  than  the  linear  treatment  that  makes  use  of  a  time  constant.  Ducoffe  and  White  (see  Bibliog¬ 
raphy)  present  an  analytical  and  experimental  study  of  pressure  response  in  a  sensing  system  which  has  a  long 
tube  and  a  reservoir  of  finite  volume.  Ducoffe 's  study  considers  typical  missile  sensing  systems  which  are  — 
subjected  to  both  impulse  and  continuous  transient  input  pressure  functions.  These  two  studies  are  concerned 
with  the  continuum  flow  regime.  Pressure  response  at  extremely  low  pressure  where  slip  and  molecular  flows 
occur  is  discussed  in  Davis  (see  Bibliography). 

The  analysis  presented  here  concerns  the  pressure  response  in  a  sensing  system  which  consists  of  a  long 
tube  connected  to  a  negligible -volume  reservoir.  This  system  may  be  represented  simply  by  a  tube  or  duct 
which  is  closed  at  one  of  its  extremities.  An  equation  is  derived  that  approximately  describes  pressure  re¬ 
sponse  in  laminar  flows  of  a  compressible  fluid  in  a  constant -diameter  tube.  The  analysis  makes  use  of  the 
Poiseuille  equation  for  laminar  flow  with  isothermal  state  changes.  A  common  derivation  of  this  equation  is 
included,  and  the  assumptions  and  limitations  in  the  analysis  are  discussed  with  reference  to  the  conditions  en¬ 
countered  during  flight  through  the  earth's  atmosphere.  Numerical  solutions  of  the  equation  are  obtained  for 
three  types  of  boundary  conditions:  step  functions,  ramp  (linear)  functions,  and  the  functions  that  represent 
constant -velocity  flight --both  ascent  and  descent --in  an  isothermal  atmosphere. 

The  first  two  types  of  functions  may  be  used  with  systems  subjected  to  impulse  pressure  functions  such  as 
blast  waves  or  shock  waves.  The  last  type  may  be  used  to  determine  pressure  response  during  variable -velocity 
flight  in  real  atmospheres. 


Analysis 

The  equation  of  continuity  may  be  derived  from  the  law  of  conservation  of  mass  using  the  concept  of  a 
continuum.  This  concept  is  applicable  when  the  me  an- free  path  X  of  the  fluid  molecules  is  small  compared  with 
the  smallest  significant  dimension  in  the  problem;  i.  e.  ,  the  tube  diameter  d  in  the  present  study.  Since  the 
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mean-free  path  in  the  earth's  atmosphere  is  less  than  2  x  inches  at  altitudes  up  to  10^  feet,  the  continuity 
equation  based  on  the  concept  of  a  continuous  fluid  is  valid  up  to  high  altitudes. for  measuring  systems  of  gen¬ 
eral  concern.  In  cylindrical  coordinates,  the  equation  of  continuity  of  a  compressible  fluid  is 


dp 

at 


r" 


ar 
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Bpvg  apv^. 
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ax 


=  0  , 


(1) 


where 

p  is  fluid  density 
t  is  time,  and 

Vq,  and  are  the  radial,  tangential,  and  axial  velocity  components  with  respect  to  the  tube  along  the 
directions  r,  0,  and  x,  respectively. 

The  equation  of  state  of  a  gas  is 


ZRT 


(2) 


where,  in  general,  the  dimensionless  compressibility  factor  Z  =  Z(p,  t).  By  substituting  Equation  2  into  Equa¬ 
tion  1  and  by  assuming  a  perfect  gas  (Z  =  constant),  the  following  equation  is  obtained 


at 


p 


+ 


(3) 


With  air  at  sea-level  pressure,  Z  is  constant  within  1/2  percent  over  the  temperature  interval  -100®  <  T 

-2 

<  4400®F,  and,  at  a  pressure  of  10  atmospheres,  Z  varies  about  1/2  percent  over  the  interval  -300®  <  T  < 
3600®F,  Thus,  over  a  broad  range  of  thermodynamic  conditions,  no  large  variations  occur  in  the  compressi¬ 
bility  factor. 

The  equations  of  motion  (Navier- Stokes  equations)  are  derived  from  Newton's  second  law  and  Stokes  law 
of  friction.  In  the  derivation  of  Stokes  law,  the  assumption  is  made  that  fluid  stresses  are  proportional  to  the 
rate  of  change  of  strain  with  time.  This  assumption  is  valid  for  laminar  flow  only.  Although  laminar  flow  has 
been  maintained  to  Reynolds  numbers  as  large  as  20,  000  by  minimizing  flow  disturbance  at  the  tube  inlet,  ex¬ 
periments  have  shown  that  steady  flow  in  a  tube  is  normally  laminar  at  Reynolds  numbers  up  to  approximately 
2000.  Further,  experiments  have  shown  that  the  motion  is  stable  in  laminar  flow--i.  e. ,  an  initially  turbulent 
flow  will  not  be  maintained  indefinitely.  The  three  equations  of  motion  of  a  compressible,  viscous  fluid  are 
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In  general,  the  viscosity  is  expressed  by  an  empirical  relation.  The  terms  F^,  F^,  and  F^  are  the  r,  0,  and 
X  components  of  the  body  forces  per  unit  volume.  In  the  absence  of  body  forces  other  than  those  due  to  gravity, 
the  magnitude  of  the  resultant  of  the  three  forces  is  F  =  pg  ,  where  g  is  the  apparent  acceleration  of  gravity. 
When  the  coordinate  system  is  in  motion,  the  apparent  gravity  may  be  much  larger  than  the  standard  gravity,  g. 
Consequently,  the  induced  pressure  gradient  can  be  significantly  large.  With  the  above  equations  it  can  be 
shown  that  when  the  fluid  is  acted  upon  by  the  body  force  only,  the  ratio  of  the  pressures  at  the  two  ends  of  a 
tube  of  length  L  is  approximately 


P2  ^  /  -1-1 

—  ^  exp  (g  R  T  L)  sin 
Pi  ^  V^x 


where 

ij/  is  the  angle  between  the  longitudinal  axis  of  the  tube  and  a  horizontal  plane, 

T  is  temperature, 

2  2 

R  is  the  gas  constant  (for  air,  R  =  1716  ft  /sec  /^Rankine),  and 
g  is  the  apparent  (absolute)  acceleration  in  the  direction  of  the  tube  axis. 

The  effect  of  the  body  force  can  be  large  when  the  tube  axis  is  parallel  to  the  direction  of  motion.  However, 

2 

even  when  acceleration  in  that  direction  is  30  times  the  standard  gravity,  i.  e.  ,  when  g^  =  30(32.  174  ft/sec  ), 
the  pressures  at  the  ends  of  a  10-foot  tube  differ  by  only  1  percent  if  T  =  lOO^F.  At  higher  temperatures,  the 
differential  pressure  is  correspondingly  smaller.  Thus,  the  effect  of  the  body  force,  F,  is  normally  small  and 
will  be  neglected  in  Equations  4,  5,  and  6. 

There  are  now  four  equations  (3  through  6)  and  five  unknown  quantities  (p,  T,  v^,  v^,  and  v^).  The  energy 
equation  is  the  additional  relation  which  is  required  to  define  the  motion.  However,  if  the  flow  process  is  iso¬ 
thermal,  there  are  only  four  unknowns  and  the  above  equations  are  sufficient.  Unfortunately,  even  with  this 
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simplification  there  is  no  known  method  for  general  solution.  The  solution  is  apparently  possible  with  the  aid 
of  a  computer,  but  further  simplifications  will  be  made  to  provide  equations  which  are  more  manageable. 

If  V  =  V  =  8v  /9x  =  9v  /90  =  0,  which  implies  a  steady,  fully  developed,  incompressible  flow  in  a 
r  y  X  X 

straight  tube  of  constant  diameter,  then  only  the  last  of  the  three  motion  equations  exists.  It  reduces  to  the 
following  expression  when  the  fluid  viscosity  is  constant  and  the  body  force  is  zero 


dx 


\r/  dr 


(7) 


For  the  known  boundary  condition  “  0  when  r  =  r^,  which  requires  that  no  slip  flow  occur  at  the  walls  and 
therefore  implies  a  continuum,  the  solution  to  Equation  7  is 


V 


X 


(8) 


where  the  axial  pressure  gradient  dp/dx  e  (p^  -  p^)/L,  and  p^  and  p^  are  constant  pressures  at  either  of  the 
open  ends  of  a  tube  of  length  L.  The  average  velocity  v  is 


V 


'  2\ 

i8ju/dx  ‘ 


(9) 


This  is  the  well-known  Hagen- Poiseuille  equation.  A  compressible  flow  will  be  approximated  by  applying  the 
Hagen- Poiseuille  expression  to  each  increment  of  tube  length.  The  average  velocity  becomes 


V  = 


-2.)ip 

\8n/  dx  ' 


(10) 


and  the  average  velocity  gradient  along  the  tube  axis  is 


3x=  (^o)a^p 

ax  \8M/g^2  • 


(11) 


By  assuming  a  quasi- steady  flow  and  substituting  Equations  10  and  11  into  Equation  3,  the  following  relation  is 
obtained 


(Tj'ar^^ra?'  1"T‘;‘^+P"a7+  iT/'ae 


© 


(g) 


m  9e3I+  p 

\T/  8x  8x  ^ 


(12) 
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Consistent  with  the  initial  assumptions  of  an  isothermal,  fully  developed  flow,  the  terms  v^,  8T/at,  3T/ax,  and 
all  derivatives  with  respect  to  r  and  0  are  zero.  Equation  12  then  reduces  to  the  following  second-order,  non¬ 
linear  equation 


at  \64)u/ 


(13) 


where  d  is  the  tube  diameter. 

In  general,  the  boundary  and  initial  conditions  for  a  tube  which  has  the  closed  end  at  x  =  L  are 

(i)  p(0,  t)  =  f(t), 

(ii)  p(x,  0)  =  f(0),  and 

(iii)  |E(L,  t)=0. 

The  first  boundary  condition  is  the  forcing  or  driving  function  and  stPtes  that  the  pressure  at  the  open  end 
(x  =  0)  of  the  tube  is  a  specified  function  of  time.  This  may  be  the  unknown  function  in  some  problems  and,  in 
that  event.  Equation  i  may  be  replaced  by  the  known  function  p(L,  t)  =  g(t).  Equation  ii  states  that  the  pressure 
is  initially  constant  throughout  the  tube.  The  last  boundary  condition  is  derived  from  Equation  10  by  specifying 
that  velocity  at  the  closed  end  of  the  tube  is  zero.  This,  of  course,  implies  negligible  volume  in  the  reservoir. 

Equation  13  may  be  made  dimensionless  by  defining  a  pressure  coefficient  P  =  p/p^.  a  length  X  =  x/L, 
and  a  time  r  =  t/t^S  where  p^  is  a  constant  reference  pressure,  L  is  the  tube  length,  and  t-  is  a  positive  incre¬ 
ment  of  time.  The  final  equation,  with  a  change  in  notation,  is 

P  =  KP^  ,  (14) 

T  XX  ' 

2  2 

where  the  dimensionless  coefficient  K  =  (p^t*d  )/(64/uL  ).  The  corresponding  boundary  and  initial  conditions 
for  the  response  equation  (Equation  14)  are 

(iv)  P(0,  t)  =  F(r), 

(v)  P(X,  0)  =  F(0),  and 

(vi)  P^(l,  r)  =  0, 

where  the  dimensionless  functions  are  F(r)  =  p  ^f(t),  F(0)  =  p  ^f(O),  and  P„(l,  r)  =  Lp  P^(L,  t). 

o  o  -A  ox 


Results 

Equation  14  is  similar  in  form  to  the  nonsteady,  one -dimensional  diffusion  equation.  However,  unlike  the 
diffusion  equation,  the  dependent  variable  is  squared  on  the  right-hand  side  of  the  pressure-response  equation. 
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Although  no  general  analytical  solution  to  the  response  equation  has  been  found,  numerical  solutions  were  ob¬ 
tained  with  an  analog  computer.  The  solutions  relate  the  pressure  coefficients  at  the  two  ends  of  the  tube  in  the 
form  of  a  ratio  ^2^^!  where  the  subscripts  1  and  2  refer  to  the  open  (X  =  0)  and  closed  (X  =  1)  ends,  respec¬ 
tively.  The  ratio  of  the  pressure  coefficients  reduces  to  the  pressure  ratio  across  the  tube  since  the  constant 
reference  pressures  p^  in  the  denominators  of  the  coefficients  cancel  each  other.  Thus,  ^2^^!  ""  ^2^?!  where 
is  the  pressure  at  the  closed  end  of  the  tube  and  p^  is  the  pressure  at  the  orifice  (open  end). 

Three  forms  of  the  general  forcing  function  F(r)  s  which  is  used  in  the  boundary  and  initial  conditions 
(Equations  iv  and  v),  were  considered:  step  functions,  ramp  (linear)  functions,  and  the  conditions  that  cor¬ 
respond  to  constant-velocity  flight  in  an  isothermal  atmosphere.  The  calculations  were  made  for  conditions 
that  represent  gas  flowing  both  out  of  and  into  the  tube  with  the  pressure  ratio  varying  from  0.  5  <  P2/P2^ 
approximately. 


For  a  step  function,  the  forcing  function  becomes  F(r)  =  1  when  r  >0  and  F(r)  =  when  r<  0.  The  con¬ 
stant  pressure  coefficient  is  =  P^/p^  where  p^  and  p^  are  respectively  the  initial  (r<  0)  and  final  (r  >  0) 
pressures  at  the  tube  inlet  (X  =  0).  Figure  1  summarizes  the  results  of  the  calculations.  Ten  step  sizes  were 

taken  over  the  pressure  ratio  range  0.  5  <  p^/p.  <  1.  5.  The  pressure  ratio  is  plotted  as  a  function  of  the  di- 

2  ^2 

mensionless  time  parameter  Kr  =  (p^d  t)/64A4L  )  where  the  reference  pressure  p^  is,  as  stated  above,  the  final 
pressure  in  the  tube,  d  is  the  tube  diameter,  t  is  time,  /u  is  the  viscosity  coefficient,  and  L  is  the  length  of  the 
tube.  The  results  reveal  that  the  time  required  for  the  pressures  at  the  tube  extremities  to  become  essentially 
equal  is  nearly  independent  of  the  initial  differential  pressure,  i.  e.  ,  the  size  of  the  step.  Further,  for  a  given 
initial  differential  pressure  and  a  given  final  pressure,  the  direction  of  flow  in  the  tube  does  not  measurably  af¬ 
fect  the  time  history  of  differential  pressure  (p^  -  p^). 


For  a  ramp  function,  the  forcing  function  is  F(r)  =  =  1  +  Ar  where  A  is  an  arbitrary  constant.  The 

pressure  coefficient  is  =  Pi^Pq  where  Pj,  =  p^[l  +  (A/t*)t]  and  p^  is  now  defined  as  the  initial  (r  =  0)  pressure 
in  the  tube.  Figure  2  shows  the  results  of  the  calculations  when  the  pressure  decreases  with  time.  For  con- 
venience,  the  value  of  A  was  selected  as  A  =  -5.  The  pressure  ratio  across  the  tube  is  plotted  versus  dimen- 

O 

sionless  time  r=  t/t^,  and  various  curves  which  represent  different  values  of  the  parameter  K  =  (p  t*d  )/ 

2  ® 

(64^L,  )  are  shown.  With  any  initial  pressure  p^,  the  pressure  ratio  across  the  tube  continually  increases  as 

time  increases.  Figure  3  is  a  plot  of  the  same  parameters  used  in  Figure  2,  but  corresponds  to  the  case  where 

the  pressure  at  the  open  end  of  the  tube  increases;  p^  =  p^[^l  +  (A/t*)tJ  where,  in  the  calculations,  A  =  50. 

The  results  show  that  with  any  initial  pressure,  the  pressure  ratio  across  the  tube  first  decreases  and  later  in¬ 


creases. 


In  an  isothermal  atmosphere,  the  ambient  pressure  p^  at  altitude  h  is  p^  =  p^^  exp  (-gR  ^T^^h)  where 
g  is  the  gravity  constant,  R  is  the  gas  constant,  T  is  the  atmospheric  temperature,  and  p„  is  the  atmospheric 
pressure  at  sea  level.  When  a  vehicle  is  ascending  in  an  isothermal  atmosphere  at  a  constant  velocity  V,  the 
pressure  p  at  any  time  t  and  any  point  on  the  vehicle  is  approximately  p  =  np^,^  exp  (-gR~^T"^  Vt  sin  0)  where 
n  is  a  known  constant  and  0  is  the  flight  path  angle  with  respect  to  a  horizontal  plane.  If  the  open  end  of  the  tube 

is  connected  to  an  orifice  on  the  vehicle  surface,  then  p  =  p^  and  the  pressure  coefficient  P.,  becomes  P.  = 
r-l-ln  ^  ^  ^ 

Pj^Pq  ~  ®^P  h*)j  where  the  reference  altitude  h*  =  Vt*  sin  0.  For  convenience  the  reference  pres¬ 

sure  is  taken  to  be  the  initial  (t  =  0  and  therefore  h  =  0)  pressure  in  the  tube;  thus  p  =  np„,  .  Hence,  the  forc- 

f  1 

mg  function  becomes  F{r)  =  Pj  =  exp  (-Cr)  where  the  constant  coefficient  is  C  =  gR  Vt*  sin  0.  Similarly, 
it  can  be  shown  that  when  the  vehicle  is  descending  at  a  constant  velocity  from  an  altitude  h*,  the  forcing 
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Figure  3.  Pressure  Response  to  Linear  Forcing  Functions  (Increasing  Pressure) 
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function  is  F{r)  =  =  exp  (Cr)  where,  again,  the  reference  pressure  is  the  initial  (t  =  0  and  h  =  h*)  pressure  in 

the  tube;  i.  e. ,  For  both  the  ascending  and  descending  flight  conditions,  the  constant  C  has 

been  taken  as  10  in  the  calculations.  At  normal  pressures  and  temperatures,  this  value  corresponds  to  an  alti¬ 
tude  range  from  sea  level  to  a  maximum  of  several  hundred- thousand  feet. 

Figure  4  shows  the  solution  to  Equation  14  with  the  boundary  condition  F(r)  =  exp  (-Cr)  and  the  initial  con¬ 
dition  F(0)  =  1.  The  pressure  ratio  p^/p  is  plotted  as  a  function  of  a  dimensionless  time  r  +  r  =  (t/t*)  +  r  , 
and  various  curves  are  shown  which  represent  arbitrary  values  of  the  parameter  K  =  (p^t^d  )/64azL  ).  A  sliding 
time  scale  is  used  with  each  curve  since,  after  an  initial  transient  phase,  an  approximate  steady-state  solution 
exists  which  is  common  to  all  values  of  K.  The  origin  for  each  of  the  response  curves  is  displaced  by  an 
amount  r  which  depends  only  on  the  value  of  K,  and  is  indicated  at  several  places  in  Figure  4.  The  value  of  r 
IS  approximately  defined  by  the  relation  10^^  ~  -ln(K/10  ).  In  the  ascending  flight  condition,  pressure  ratio 
first  increases  at  a  rapid  rate,  and,  after  the  transient  phase,  the  rate  of  increase  becomes  relatively  small 
(K  »  1).  However,  the  pressure  ratio  continues  to  increase  as  altitude  (time)  increases.  The  solutions  to  the 
response  equation  indicate  that  even  if  the  tube  is  short  and  the  flight  velocity  is  low,  eventually  an  altitude  will 
be  reached  above  which  the  pressure  ratio  will  exceed  an  arbitrary  value.  Of  course,  the  mean-free  path  of 
the  gas  in  the  tube  normally  decreases  with  altitude,  and  when  that  dimension  becomes  comparable  with  the  tube 
diameter,  the  continuum  concept  is  no  longer  satisfactory.  The  criterion  for  this  limitation  is  considered  later 
in  the  "Discussion”  section. 

o 

Note  that  the  application  of  the  results  is  not  limited  to  values  of  the  parameter  K  less  than  10  which  is 

3 

the  maximum  value  shown  in  Figure  4.  If  K>  10  ,  then  <  0,  and  the  pressure  ratio  is  approximately  unity 
during  the  time  interval  0  <  r  <  I  ^^1  •  In  this  time  interval,  the  pressure  ratio  or  differential  pressure  cannot 
be  determined  without  cross  plotting  and  extrapolating  the  calculations.  When  r>  I  the  existing  data  are 
sufficient  to  determine  pressure  response. 

Figure  5  shows  the  solution  to  Equation  14  with  the  descending- flight  boundary  conditions;  F(r)  =  exp  (Cr) 
where  C  =  10,  and  F(0)  =1.  As  in  the  previous  figure,  a  sliding  time  scale  is  used  with  each  curve.  The  pres¬ 
sure  ratio  is  plotted  versus  the  dimensionless  time  ^  where  is  approximately  defined  by  the  relation 
lO^o  ^  ln(2K).  The  results  are  similar  to  those  corresponding  to  the  ascending  flight  conditions  in  that,  after 
the  transient  phase,  an  approximate  steady- state  solution  exists  which  is  common  to  all  values  of  the  parameter 
K.  Thus,  these  data  may  be  used  to  determine  the  pressure  response  of  a  measuring  system  in  a  vehicle  which 
is  within  the  applicable  altitude  range,  even  though  the  vehicle  may  have  descended  from  higher  altitudes  where 
the  continuum  concept  is  not  necessarily  valid  or  where  the  initial  condition  differs  from  that  used  in  the  calcu¬ 
lations. 

The  results  that  are  shown  in  Figures  4  and  5  may  be  used  to  determine  the  approximate  pressure  re¬ 
sponse  of  sensing  systems  in  real  atmospheres.  This  requires  iteration  since  the  constant  K  is  variable  in  a 
real  atmosphere.  The  reference  pressure  p^,  which  is  used  in  the  calculation  of  K,  will  be  a  fictitious  initial 
pressure  calculated  for  an  isothermal  atmosphere  of  temperature  T^,  but  reduces  to  a  true  pressure  in  the 
real  atmosphere  at  the  desired  altitude.  If  the  vehicle  is  accelerating  or  decelerating  as  will  be  the  case  with 
general  trajectories,  the  pressure  coefficient  will  be  a  variable  and  n  will  not  be  constant.  This  effect  may 
be  taken  into  account  by  dividing  the  flight  path  into  parts  that  are  approximately  constant -velocity  intervals, 
and  applying  the  calculating  procedure  to  each  interval. 
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Figure  4.  Pressure  Response  During  Ascending  Flight  in  Isothermal  Atmospheres 
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Discussion 


As  indicated  previously,  the  analysis  is  applicable  only  when  the  Knudsen  number  is  much  less  than  unity. 

Thus,  there  is  some  small  Reynolds  number  below  which  the  the  response  equation  is  not  valid.  The  mean-free 

2  1/2 
path  A  can  be  written  A  =  kT  /(^p^ttu  ),  where  k  is  Boltzmann's  constant,  T-  is  temperature,  ^  =  2  '  if  the 

11  2 

velocity  distribution  is  Maxwellian,  and  the  product  ira  is  the  collision  cross-sectional  area  for  a  rigid  spher¬ 
ical  molecule.  If  T^  is  measured  in  degrees  Rankine  and  p^  is  pounds  per  square  foot,  the  mean-free  path  of 

air  becomes  A  =  8.  88  x  10"'^(T  /p. )  feet.  Since  the  dimensionless  Knudsen  number,  A/d,  must  be  small  com- 

1-7  ^ 

pared  with  unity,  then  8.  88  x  10  (T  /p.d)  «  1.  Therefore,  the  response  equation  is  valid  for  flight  conditions 

11  _  5 

in  the  earth's  atmosphere  when  p^d/T^  »  10  ,  where  p^  is  pressure  in  pounds  per  square  foot,  d  is  tube  di¬ 

ameter  in  inches,  and  T^  is  temperature  in  degrees  Rankine. 

In  actual  flows,  the  radial  distribution  of  velocity  is  not  fully  developed  near  the  tube  inlet.  The  length  ^ 
required  to  attain  fully  developed,  laminar  flow  is  a  function  of  Reynolds  number  and  is  discussed  in  Schlicting 
(see  Bibliography);  the  required  length  in  a  steady,  laminar  flow  is  approximately  i/d  =  0,  03  R^^,  where 
Reynolds  number  is  based  on  tube  diameter.  Since  in  the  analysis  presented  here  a  fully  developed  profile  was 
assumed,  it  is  implied  that  the  tube  length  L  is  large  compared  with  the  dimension  i.  Therefore,  the  response 
equation  is  limited  to  the  geometry  range  L/d  »  0.  03  R^. 

Shapiro  (see  Bibliography)  shows  that  there  is  a  limiting  Mach  number  M  for  continuous  isothermal  flow 

-1/2 

in  a  tube.  The  magnitude  of  the  limiting  Mach  number  in  a  steady  flow  is  M  =  7  ,  where  7  is  the  ratio  of 

specific  heats.  When  a  subsonic,  isothermal  flow  approaches  the  limiting  condition,  choking  effects  occur  and 
fluid  properties  vary  rapidly.  It  is  expected  that  the  response  equation  will  not  describe  the  pressure  distribu¬ 
tion  accurately  near  the  limiting  condition,  since,  in  this  analysis,  the  possibility  of  large  variations  in  fluid 
properties  was  precluded.  Further,  Shapiro  shows  that  to  maintain  isothermal  flow  at  the  limiting  condition, 

an  infinite  he  at -transfer  rate  is  required.  Hence,  the  basic  assumption  of  isothermal  flow  is  not  realistic  when 
-1/2 

M  -►  7  .  The  application  of  the  results  presented  in  Figures  1  through  5  should  be  therefore  limited  to  the 

corresponding  pressure -ratio  range. 

No  experiments  have  been  made  and  there  are  no  available  data  to  compare  with  the  analytical  results. 
However,  when  the  flow  is  steady,  the  response  equation  can  be  integrated  to  define  the  axial  pressure  distribu¬ 
tion,  and  a  comparison  can  be  made  with  existing  experimental  data.  This  comparison  will  provide  an  indica¬ 
tion  of  the  accuracy  of  the  response  equation. 

If  the  flow  is  steady,  then  =  0,  and  Equation  14  reduces  to  the  following  differential  equation 


(15) 


With  the  two  boundary  conditions 

(vii)  P(0)  =  P^  , 

(viii)  P(l)  =  P2, 
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Equation  15  can  be  easily  integrated.  The  result  is 


P/P^ 


(P2/P1) 


(16) 


where  P  -  P(X),  X  =  x/L.  Since  the  pressure  coefficient  is  based  on  a  constant  reference  pres¬ 

sure^  Equation  16  reduces  to 


P/P^  = 


1  + 


(P2/P1) 


[x/l] 


1/2 


(17) 


where  and  are  respectively  the  inlet  (x  =  0)  and  exit  (x  =  L)  pressures.  This  relation  describes  the  axial, 
static  pressure  distribution  in  a  constant-diameter  tube  and  it  pertains  to  laminar,  isothermal  flow  of  a  com- 
pressible  fluid. 

Equation  17  can  be  also  derived  from  steady-flow  considerations.  That  type  of  derivation  is  reported  in 
Ducoffe  s  study  (see  Bibliography)  and  the  equation  is  compared  with  the  results  of  experiments  using  air.  In 
the  experiments,  exit  pressure  was  approximately  equal  to  ambient  pressure  at  sea  level.  Consequently,  the 
Reynolds  number  is  large;  the  maximum  Reynolds  number,  based  on  tube  diameter,  is  of  order  10^.  The 
pressure  ratio  p/p^  is  plotted  as  a  function  of  axial  station  in  Figures  6  and  7,  and  a  comparison  is  made  of  the 
experimental  data  with  Equation  17.  In  Figure  6,  which  pertains  to  a  tube  with  a  length-diameter  ratio  of  667, 
the  experimental  data  and  the  theoretical  estimate  differ  by  a  few  percent.  In  Figure  7,  which  corresponds  to  a 
length-diameter  ratio  of  229,  the  agreement  is  within  a  few  percent  only  when  the  pressure  ratio  is  approxi¬ 
mately  unity.  Equation  17  and  experimental  data  are  in  best  agreement  when  the  length-diameter  ratio  of  the 
tube  is  large  and  the  Reynolds  number  is  small  (R^  <  2000). 


Conclusion 

An  equation  is  derived  that  approximately  describes  pressure  response  in  a  system  that  consists  of  a  long 
tube  connected  to  a  reservoir  of  negligible  volume.  The  equation  pertains  to  the  laminar  flow  of  a  compressible 
fluid  in  a  constant -diameter  tube  and  is  applicable  in  the  continuum  regime.  Numerical  solutions  of  the  equa¬ 
tion  are  presented  for  three  types  of  boundary  conditions:  step  functions,  ramp  functions,  and  the  functions  that 
represent  ascent  and  descent  flight  at  constant  velocity  in  an  isothermal  atmosphere.  These  solutions  may  be 
used  to  estimate  the  pressure  ratio  or  differential  pressure  across  the  tube  during  transient  conditions  of  flight 
in  real  atmospheres.  They  may  also  be  used  to  determine  the  response  of  systems  subjected  to  impulse  pres¬ 
sure  functions. 
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